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ABSTRACT 

We consider the Renormalization-Group coupled equations for the effective 
potential ^(</>) and the field strength Z^cf)) in the spontaneously broken phase as 
a function of the infrared cutoff momentum k. In the A; ^ limit, the numerical 
solution of the coupled equations, while consistent with the expected convexity 
property of ^(0), indicates a sharp peaking of ^(0) close to the end points of 
the flatness region that define the physical realization of the broken phase. This 
might represent further evidence in favor of the non-trivial vacuum field renor- 
malization effect already discovered with variational methods. 



Pacs ll.lO.Hi , ll.lO.Kk 



1. Renormalization Group (RG) techniques originally inspired to the Kadanoff- 
Wilson blocking procedure P represent a powerful method to approach non- 
perturbative phenomena in quantum field theory. A widely accepted technique 
consists in starting from a bare action defined at some ultraviolet cutoff A and 
effectively integrating out shells of quantum modes down to an infrared cutoff k. 
This procedure provides a fc— dependent effective action Tk[^] that evolves into 
the full effective action r[$] in the limit 0. 

The fc— dependence of ^k[^] is determined by a differential functional flow 
equation that is known in the literature in slightly different forms |2l El El El E] ■ 
In particular, with the flows discussed in detail in Ref.|,.7^ one starts form flrst 
principles and obtains a class of functionals that interpolates between the classical 
bare Euclidean action and the full effective action of the theory. However, some 
features, such as the basic convexity property of the effective action for k ^ 
[HI Ellini IHj , are independent of the particular approach. 

In this Letter, we shall study the coupled equations for the fc— dependent 
effective potential Vk{4>) and fleld strength Zk{(f)), which naturally appear in a 
derivative expansion of ^k[^] (around the coordinate independent conflguration 
$(x) = ), using a proper-time infrared regulator. This approach generates 
evolution equations which are well deflned truncations of the flrst-principle flows 
within a background fleld formulation ^21 El- Incidentally, as shown in they 
are also an approxmation to another exact flow (a generalised Callan-Symanzik 
flow) without background flelds. 

The equations for Vk{4>) and Zk{(f)), derived in Ref . [Tlj . correspond to such 
first-principle fiows up to additional terms explicitly determined in Ref. |12j. An 
indication of the reliability of such an approximation is provided by the consistent 
computation of the critical exponents in scalar self-interacting theories for various 
numbers D of the space-time dimensions and N of field components At the 
same time, in Ref.^Il, it was shown that going beyond the approximation Z = lis 
essential to reproduce successfully the energy gap between the exact ground state 
and the first excited state of the double well potential in the quantum-mechanical 
limit of the theory D = 1. In addition, the coupled equations obtained by using 
the proper time regulator are not affected by singularities and/or ambiguities 
that instead appear using a sharp infrared cutoff 15 . 

With these premises, assuming the possibility to neglect the additional terms 
of Ref. 1121 adopting the same notations of Refs. [TT| ITlj . we shall start our 
analysis considering the two equations 
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where we have set \^ = Vfc(0), Z = ^^(0) and used the notation V , Z',...to in- 
dicate differentiation with respect to (p. In the following, we shall first analyze 
these two equations to understand the approach to convexity and obtain informa- 
tions on Z and finally provide a possible physical interpretation of our numerical 
results. 



2. For a numerical solution of Eqs.(Q) and Q it is convenient to use di- 
mensionless variables defined as follows: t = ln(A/A;), x = k^'^^'^cp, V{t,x) = 
k~^Vk{4>) ^-iid Z{t,x) = Zk{(f)). By defining the first derivative of the effective 
potential f{x,t) = dxV(t,x), Eqs. (jll2j) become 

dl ^ {D + 2) {2-D) dl _ ^ 

dt 2 ^ 2 dx (47r)^/2 dx ^ ' 
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It is easy to show that these two coupled equations can be transformed into the 
structure (z,j=l-3) 

dt + " dx 

where the components of the vector Ui{x,t) are the unknown functions of the 
problem and where Pij, Qi and Ri can depend on x, t, Ui, 



In this way, the numerical solution has been obtained with the help of the 
NAG routines that integrate a system of non-linear parabolic partial differential 
equations in the [x, t) two-dimensional plane. The spatial discretisation is per- 
formed using a Chebyshev C° collocation method, and the method of lines is 
employed to reduce the problem to a system of ordinary differential equations. 
The routines contain two main parameters (the size of the discretisation grid 
-^points in the spatial dimension and the local accuracy A in the time integra- 
tion) that can be changed to control the stability of the solution. In our case, 
after reaching sufficiently large values of A/'points (~ 2000) and sufficiently small 
values of A (~ 10"''), the numerical results remain remarkably stable for further 
variations of these parameters. 

We shall focus on the quantum-field theoretical case D = 4 assuming standard 
boundary conditions at the cutoff scale : i) a renormalizable form for the bare, 
broken-phase potential 

^a(0) = -^^20' + (6) 
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and ii) a unit renormalization condition for the derivative term in the bare action 

Za(0) = 1 (7) 

With this choice of the classical potential the problem is manifestly invariant 
under the interchange —(f) at all values of k. 

We shall also concentrate on the weak coupling limit A = 0.1, fixing M = 1 
and A = 10. In this way, one gets a well defined hierarchy of scales where the 
infrared region corresponds to the limit A; -C M <^ A. 

Before addressing the full problem, we have considered the standard approx- 
imation of setting Z = 1. In this case, we can find a good approximation to 
the exact solution of Eq.(j2I) for large x and large t as a cubic polynomial with 
dependent coefficients 

/asy(x, t) = A{t)x'^ + B{t)x^ + C{t)x + D{t) (8) 

This type of form, motivated by RG arguments, see e.g. [^, becomes exact where 
/' is large and positive so that e^^ yielding A{t) = Aq, Bit) = Boexp{t), 
C{t) = Coexp(2t), D{t) = Z)o6xp(3t). Higher powers x", with n > 3, might 
also be inserted but they are suppressed by exponential terms e"*^""^)*. At the 
same time, / = is also a solution of Eq.Q and corresponds to a flat effective 
potential in x. As we shall see, the fc— evolution, for k ^ 0, turns out to be a 
sort of way to interpolate between / = and the asymptotic trend in Eq.®. 

We show in Fig.l the plots of \4'(0) = at various values of the infrared 

cutoff down to the smallest value k = 0.05 that can be numerically handled by 
our integration routine. This corresponds to a maximal value tmax ~ 5.3 which 
is comparable to the highest value t = 5 reported in Fig. la of Ref. 9J. In our 
problem, this seems to mark the boundary of the region that is very difficult to 
handle numerically. 

As one can see, the approach to convexity is very clean, in full agreement with 
the general theoretical arguments of Ref. [HI. It is characterized by an almost linear 
behaviour in the inner 0— region that matches with the outer, asymptotic cubic 
shape discussed above. In this sense, the limitation in the maximal t is just a 
numerical artifact of our integration method since from the physical point of view 
there should be no conceptual problem in the limit t — oo. This is supported by 
a fit to the slope S{k) of Vl{(j)) at = in Fig.l, which suggests the following 
functional form (a and b being numerical coefficients obtained from the fit) 

S{k) ~ a[exp{-bk^) - 1] (9) 

that vanishes in the limit k 0. 

For small k, the minimum of the first derivative of the effective potential 
(f> = (j){k) does not correspond to an analytic behaviour. This is a well known 
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result of quantum field theory [T7j: the Legendre-transformed effective potential 
is not an infinitely differentiable function. 

Let us now consider the full problem defined by Eqs.© and (jH). Again, for 
large x and t, the pair / = fasy{x,t) and Z = 1 provide a simultaneous solution. 
However, for finite values of k, where the potential is not yet convex downward 

and /' < 0, the term e z drives Z to grow. 

We show in Fig. 2, the simultaneous solutions for Zk{(j)) (upper panel) and 
= (lower panel) for various values of the infrared cutoff down to 

k = 0.07 (again tmax ^ 5) that represents, for the coupled problem, the point 
beyond which the integration routines no longer work. The strong peaking in 
Zk{4>), to a very high accuracy, occurs at 0(fc) where 14'(0) minimum. On 

the basis of the general convexification property this point tends, for A; ^ 0, to 
the end point 0(0) = 0o of the flatness region that defines the physical realization 
of the broken phase. Finally, we report in Fig.3 the values of Zk{(j)) at = 
(circles) and at the peak for = (j){k) (diamonds) in the range 0.07 < k < 0.15. 
Using some extrapolation forms, the value of the peak seems to tend, in the limit 
A; — * 0, to a very large but finite value Zk=o{(po)- 

3. Let us now explore a possible physical interpretation of the numerical re- 
sults reported above. We start by observing that spontaneous symmetry breaking 
is usually considered a semi-classical phenomenon, i.e. based on a classical po- 
tential with perturbative quantum corrections. These, with our choice of the 
bare parameters A = 0.1 and M = 1 in Eq.®, and our cutoff value A = 10, 
are typically small for all quantities. In particular Z, in perturbation theory, is a 
non-leading quantity since its one-loop correction is ultraviolet finite. Therefore, 
the deviations from unity are expected to be very small. 

Re-writing the term in the standard form with Agt = 2.4, the per- 

turbative prediction is 

Zpert - 1 = O(^) ~ 10-^ (10) 

This is also consistent with the assumed exact "triviality" property of the theory 
for D = 4 jlHI that requires Z ^ 1 in the continuum limit A ^ oo. 

Now, let us compare this prediction with our Zk{(f)) in Fig. 2. For large values of 
the infrared cutoff, when the effective potential is still smooth, we find Zk{(j)) ~ 1 
for all values of 0, as expected. However, at smaller k, say k < S with S ~ 0.15, 
there are large deviations from unity in the region of where the smooth form of 
the perturbative potential evolves into the typical non-analytical behaviour of the 
exact effective potential. This leads to the observed strong peaking phenomenon 
at the point 0(A;), where V^'(0) reaches its minimum value. This point, on the 
basis of the general convexification property of V{(j)) tends, for /c — > 0, to the 
value 0(0) = 00, the end point of the flatness region that defines the physical 
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realization of the broken phase. 

If we express the full scalar field <l>(x) as 

$(x) = + /i(x) (11) 

the above results indicate that the higher frequency components of the fluctuation 
field h{x), those with 4-momentum such that 6 < \p\ < A, represent genuine 
quantum corrections for all values of the background field in agreement with 
their perturbative representation as weakly coupled massive states. 

On the other hand, the components with a 4-momentum such that \p\ < S, 
are non-perturbative for values of the background field in the range < < 
0(|j9|). In particular, the very low- frequency modes with \p\ — >■ behave non- 
perturbatively for all values of the background in the full range < < 0o- They 
can be thought as collective excitations of the scalar condensate and cannot be 
represented as standard massive states. 

The existence of a peculiar p^ ^ limit in the broken phase, for which 
one can give some general arguments 19 , finds support in the results of lattice 
simulations of the broken-symmetry phase (see Ref. |2Uj). There, differently from 
what happens in the symmetric phase, the connected scalar propagator deviates 
significantly from (the lattice version of) the massive single-particle form l/(p^ -|- 
const) for — 0. In particular, looking at Figs. 7, 8 and 9 of Ref. [21], one 
can clearly see that, approaching the continuum limit of the lattice theory, these 
deviations become more and more pronounced but also confined to a smaller and 
smaller region of momenta near = 0. 

This observation suggests that the existence of a non-perturbative infrared 
sector in a region < \p\ < 6 might not be in contradiction with the assumed 
exact "triviality" property of the theory if, in the continuum limit, the infrared 
scale 6 vanishes in units of the physical parameter m associated with the massive 
part of the spectrum. This means to establish a hierarchy of scales 6 <^ m <^ A 
such that - ^ when x 0. 

m A 

If this happens, the region < \p\ < 6 would just shrink to the zero-measure 
set Pfj_ = 0, for the continuum theory where m sets the unit mass scale, thus 
recovering the exact Lorentz covariance of the energy spectrum since the point 
Pn = forms a Lorentz- invariant subset. In this limit, the RG function Zk{(f)) 
would become a step function which is unity for all finite values of k (and 0) 
and is only singular for /c = in the range < < 0o. In this way, one is left 
with a massive, free-field theory for all non-zero values of the momentum, and 
the only remnant of the non-trivial infrared sector is the singular re-scaling of 
(the projection of the full scalar field $(x) onto p^ = 0). 

This is precisely the scenario of Refs. |22j . where for A — oo the re-scaling 
of the scalar condensate diverges logarithmically as ~ InA and the re-scaling of 
the finite-momentum modes tends to unity. The existence of such a divergent re- 
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scaling factor for the vacuum field would have potentially important phenomeno- 
logical implications for the scalar sector of the standard model and for the validity 
of the generally accepted upper bounds on the Higgs boson mass [22] . 

Of course, for a more precise comparison with the numerical results obtained 
in this Letter, one should study the value of the peak in ^^(0) at different values 
of the bare parameters to check the predicted logarithmic behaviour Zpeak ~ ln(A) 
that, at the present, is just a conjecture suggested by previous works. In turn, 
this requires to improve on the present integration routines to extend the solution 
of the RG equations towards the point k = thus reducing the arbitrariness 
associated with different extrapolation forms. 

At the same time, the existence of a peak in Zk{(l)) is an interesting feature 
that would deserve to go beyond the present numerical analysis. In this respect, 
since, according to Vk{(f>), for /c — >• 0, is only weakly dependent on the type of 
infrared regulator, we expect the peak of Z to be a physical phenomenon and not 
an artifact of the regulator. This can be checked by re- considering the problem 
from scratch at the level of the exact flow equations. 
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Figure 1: V'k{(j)) vs. 0, with Z = 1 fixed, at various values of the infrared 
cutoff k: the lowest curve is for /c = A = 10, and then, from bottom to the top, 
A; = 0.3, 0.13, 0.1, 0.08, 0.05. 
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Figure 2: Zk{(p) (upper panel) and V'k{(f)) (lower panel), vs. (p at k —0.3, 0.15, 
0.13, 0.12, 0.1, 0.07, from bottom to the top. 
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Figure 3: Zk{4>) at the origin = (circles) and at the peak (j) = 4> (diamonds), 
at different values of the infrared cutoff k. 
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